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We present the quadratic algebra of the generahzed MICZ-Kepler system 
in three-dimensional Euclidean space and its dual the four dimensional 
singular oscillator in four- dimensional Euclidean space E^. We present their 
realization in terms of a deformed oscillator algebra using the Daskaloyannis 
construction. The structure constants are in these cases function not only of 
the Hamiltonian but also of other integrals commuting with all generators 
of the quadratic algebra. We also present a new algebraic derivation of the 
energy spectrum of the MICZ-Kepler system on the three sphere using a 
quadratic algebra. These results point out also that results and explicit for- 
mula for structure functions obtained for quadratic, cubic and higher order 
polynomial algebras in context of two-dimensional superintegrable systems 
may be applied to superintegrable systems in higher dimensions with and 
without monopoles. 



1 Introduction 

The most well-known superintegrable systems are the Kepler- Coulomb sys- 
tem and the harmonic oscillator [1,2,3]. A systematic study of superintegrable 
systems with second-order integrals of motion ('quadratically superintegrable 
systems') in two-dimensional (2D) Euclidean space was begun some years ago 
[4,5]. For a review of 2D superintegrable systems we refer the reader to [6]. 

In dimensions, the symmetry algebra generated by the integrals of mo- 
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tion of the harmonic oscillator is the Lie algebra su{N), while that of the 
bounded states of the hydrogen atom is the Lie algebra so{N + 1). A non- 
linear symmetry algebra was discussed for the first time in context of the 
Kepler- Coulomb system and harmonic oscillator on the three sphere [7]. 
A quadratic algebra for these two systems and the Hartmann potential were 
obtained for fixed energy and angular momentum [8]. The integrals of the 
Hartmann system form a quadratic Hahn algebra QH{3), while a quadra- 
tic Racah algebra QR{3) describes the symmetry algebra of the harmonic 
oscillator and Kepler- Coulomb system on a 3D space of constant curvature. 

A general quadratic algebra in classical and quantum mechanics, genera- 
ted by second order integrals of motion, was studied by Daskaloyannis [9], 
who discussed its realizations in terms of deformed oscillator algebras and 
found a method to obtain the energy spectrum algebraically. We genera- 
lized this method for systems with a second and a third order integral of 
motion, generating a cubic algebra [10]. In further work all quantum super- 
integrable systems separable in Cartesian coordinates with second and third 
order integrals were studied from the point of view of polynomial algebras 
[11], and a class of quintic and seventh order algebras, their realizations in 
terms of deformed algebras, and their finite dimensional unitary representa- 
tions, were presented. Families of systems with algebras of arbitrary order 
were also considered. We also discussed the relations between superintegrable 
systems and their polynomial sjTumetry algebras, ladder operators and their 
polynomial Heisenberg algebras and supersymmetry [11]. 

A systematic search for superintegrable systems with magnetic fields was 
started in recent years [12-15]. However, the problem of degeneracy in the 
presence of a magnetic monopole was studied in earlier articles independently 
by Zwanziger [16] and by Mcintosh and Cisneros [17]. This system and his 
dynamical symmetry algebra was also studied in Ref.18. This 'MICZ-Kepler' 
system, of a magnetic monopole in fiat 3D space with an arbitrary electric 
charge plus a fixed inverse-square term in the potential, is superintegrable, 
and has an so(4) symmetry algebra generated by the Poincare vector together 
with a Runge-Lenz vector. The MICZ-Kepler problem also exists in higher 
dimensions, just as the Kepler problem does [19], where they remain su- 
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perintegrable. The Yang-Coulomb system (YCS) is another superintegrable 
system, consisting of the Yang monopole and a particle coupled to the mo- 
nopole by isospin su{2) and the Kepler- Coulomb interaction [20,21]. This 
system is related to the 8D harmonic oscillator and has a so(6) symmetry 
algebra [22]. 

The MICZ-Kepler problem was also considered in 5*^ [23]. The integrals of 
motion of this superintegrable system form a cubic algebra similar to that 
of the Kepler-Coulomb problem on S^. Further, a generalized MICZ-Kepler 
system was introduced by Mardoyan [24]. This system can be seen as the 
intrinsic Smorodinsky-Winternitz system with monopole in 3D Euclidean 
space. It can be transformed into a 4D singular harmonic oscillator by a dua- 
lity transformation [25]. The corresponding Schrodinger equation was solved 
for these two cases, but no polynomial algebras were found. 

In light of these results one can ask the question if a polynomial algebra 
can be generated by the integrals of motion of the generalized MICZ-Kepler 
system and the 4D singular harmonic oscillator. The purpose of this paper is 
to obtain such a polynomial algebra and its realization in terms of a deformed 
oscillator algebra of these two systems, and to present a new derivation of 
their energy spectrum using an algebraic method. 

The paper is structured as follows. In Section 2, we recall results on qua- 
dratic algebras and their realization in term of deformed oscillator algebras. 
In Section 3 we consider the case of the generalized MICZ-Kepler system. We 
present the quadratic algebra, the Casimir operators, the realization in terms 
of a deformed oscillator algebra and its spectrum. In Section 4, we consider 
the case of the 4D singular oscillator. In Section 5 we present a new derivation 
of the energy spectrum of the MICZ-Kepler system on by constructing 
a quadratic algebra and using its realization in terms of deformed oscillator 
algebras. 
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2 Quadratic algebra and realizations in term of 
deformed oscillator algebras 



The most general quadratic algebra [7] generated by the integrals of motion 
of a quadratically superintegrable system is 

[A,B]=C, 

[A,C] = PA^ + -i{A,B} + 5A + eB + C, (2.1) 

[S, C] = aA^ - -fB^ - B}^dA - bB^ z. 

In the [9-11], the structure constants of the polynomial algebra were assumed 
to be polynomial functions of the Hamiltonian only. Such algebras were stu- 
died in the context of 2D superintegrable systems. However, because we are 
interested in application in context of superintegrable systems in dimension 
higher than two, in this paper we allow the structure constants /3, (5, 7, e, C, 
a, d and z to be polynomial functions not only of the Hamiltonian but also 
of any other integrals of motion (Fj) that commute with the Hamiltonian, 
with each other and also with the generators of the quadratic algebra A, B 
and C (that is, {Fi,U\ = [Fi, Fj] = [F„A] = [F„ B] = [Fi,C] = 0). When 
we study this algebra's realization in terms of deformed oscillator algebras 
and its representations, we will fixed the energy [Hip = Eip) but also these 
other integrals [Fii/j = fiip) of motion that form, with the Hamiltonian, an 
Abelian subalgebra. The Casimir operator ([iT, A] = [i^, -B] = [-^, C*] = 0) of 
this quadratic algebra is thus given in terms of the generators by 

K = C^- a{A\ B} - 7{A, B^} + (^7 - b){A, B} + (7^ - e)B^ (2.2) 

+ (75 - 2C)B + + (d + ^ + a^)A^ + (^ + a6 + 2z)A. 

The Casimir operator K will be also rewritten as a polynomial of H and F,. 

The realizations of the quadratic algebra in terms of deformed oscillator 
algebras were obtained in [9] : 

[N,h^] = h\ [N,h] = -h,hh^ = ^N + l),h^h = ^N), (2.3) 
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and the structure function for the case e = 5 = (3 = 0is given by 

$(Ar) = -so72-f^K{-l+2{N+u)y-48-f%-d-f'^){-3+2{N+u)){-l+2{N+u))\l+2{N+u)) 

(2.4) 

+7^(4a7)(-3+2(iV+n))2(-l+2(iV+M))^(l+2(iV+n))2+768(472c)2+327^(-l+2(iV+n))2 

+ (-1 - 12{N + u) + 12{N + ii)2)(87^2) - 25Q-f\-l + 2(A^ + u))\-4-f^z). 

(Such reahzations for cubic and higher order polynomial algebras have also 
been obtained [10,11].) 

To obtain unitary representations we should impose the following three 
constraints on the structure function : 

^{p+l,u,k) = 0, $(0,M,A;) = 0, <j){x)>0, Vx>0 . (2.5) 

The energy E and the arbitrary constant u are solutions of the equation 
obtained from these constraints. 

3 Generalized MICZ-Kepler systems 

The generalized MICZ-Kepler system was introduced and studied by Mar- 
doyan [24], where it was shown that the Schrodinger equation for this system 
allows separation of variables in spherical and parabolic coordinates, and 
the spectral problem was solved. We use the system of units for which ( 
h = m = e = c= l) 

The Hamiltonian of the generalized MICZ-Kepler system is given by 

^=^HV-5A7 + ^-- + -^ + -^, (3.1) 
2 2r^ r r(r + z) r(r — z) 

with ^ 

A=- r(y,-x,0). (3.2) 

r(r — z) 

Let us recall the Poincare and Laplace-Runge-Lenz vectors for the MICZ- 
Kepler system on the three-dimensional Euclidean space E3 : 

J = fx (-zV - sA) - S-, (3.3) 
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/' = 1(J X (-zV - sA) - (-zV - sA) X /+ -. (3.4) 
2 r 

The generalized MICZ-Kepler system has the following integrals of motion 

{IH,A] = [H.B] = \H,J,]=0) : 

A=P + j!^ + Jlfi^, B = + _ £f±4, (3.5) 

[r + z) (r — 2;) r(r + 2;) r(r — z) 

From these, we can form the quadratic algebra (with [A, J^] = [B, Jz] = 0) 

[A,B] = C, 

[A, C] = 2{A, B} + AsJz + 4(c2 - ci), (3.6) 

[B, C] = -2B^ + 8HA - UlH + 4(1 - - 2ci - 2c2)H + 2, 

which is the Daskaloyannis quadratic algebra at the special values of its 
structure constants 

/3 = 5 = e = a = 0, 7 = 2, d = m, C = AsJz + A{d-c), (3.7) 

z = {-AJl - As^ - 4(2ci + 2c2 - l))H + 2. 

The Casimir operator can be rewritten in terms of the Hamiltonian and 
(which commute with all generators of the quadratic algebra) as 

K= -Ss^JlH + lQ{c-d)sJzH-S{ci-C2fH + AJl + A{2ci + 2c2 + s^). (3.8) 

Using (2.4), (3.7) and (3.8) (with the energy and the z component of the 
Poincare vector taking the values Hip = Eip and Jzip = mil)) we then obtain 
the structure function (with mi = ^J{m — sY + 4ci, m2 = y^^?T7r+"sp~T4c2) 

$(x) = -3 ■ 2''^Eix + u-i- + —^)){x + n - (- - ^^)) (3.9) 

1 .mi + m2..., A .mi-m2 
(x + M - (- - ( ^ ))){x + n - (- - ( ^ ))) 

1 ,-mi + m2.-.. ,1 .-mi-m2,-, 
(x + w - (- - ( ^ )))(x + w - (- - ( ^ ))) 
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From the constraint (2.5) we then obtain the following finite dimensional 
unitary representation and energy spectrum 

E=— — -, u = - + ^^=, (3.10) 

2(p+ 1 + ?Tii + m2)2' 2 7^2^' ^ ' 

^{x) = 3 ■ 2' \ <P+^-^) (p + 1 + _ + 1 + ^2 - x) (3.11) 
(p + 1 + mi + m2)^ 

{p + 1 + nil + ^2 ~ {2p + 2 + nil + '>t^2 — x) . 

In the Ref.23 the energy spectrum in terms of the principal quantum number 
was presented 

E = Ei'^ = ^i-^ — (3.12) 

2(n+^)2 ^ ' 

The Eq.(3.9) of the Ref.24 relate the parabolic quantum numbers ni and n2 
with the principal quantum number n 

\m — s\ + \m + s\ ^ , , ^ , , J. 

n = ni+n2-\ hi, nii = \ni — s\ + Oi, ni2 = \ni + s\ + 02 

(3.13) 

By taking p = ni + n2 the formula given by Eq.(3.10) and Eq.(3.12) coincide. 

To summarize : We have obtained the quadratic algebra for the generalized 
MICZ-Kepler system, constructed its realization in term of deformed oscilla- 
tor algebra, and obtained the energy spectrum. This is the first study of this 
system from the point of view of its polynomial algebra. 



4 Four dimensional singular oscillator 

A relation between systems with and without monopoles has been discussed 
in [22], and such a relation was obtained for the generalized MICZ-Kepler 
system system and a 4D singular harmonic oscillator [25,26], which is also 
superintegrable. Let us recall results concerning these systems and the duality 
transformation. We consider the following change of variables [25] : 

tp^'\r) f-^ ^(f,7) = tp'^'\r) ^ , , s^-i—, 7e[0,47r). (4.1) 

V47r (77 
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The generalized MICZ-Kepler system is then transformed into a 4D double 
singular oscillator [26]. If we change /3 = 9 and a = (p and pass from the 
coordinates r,a, /3,^ to the coordinates 

^ ■ Q!-|-'y ■ ct — 'y 

Mo + mi = ncos( — )e*^~, U2 + iu^ = usm{—)e^^~ , (4.2) 
we obtain the Hamiltonian 

^ = -o(7r^ + 7r^ + 7r^ + 7r^)+^ + 7^^ + 7^^- (4-3) 



This system was studied in [26], and is separable in Eulerian, double-polar 
and spheroidal coordinates in E4. It has the following integrals of motion ( 

[H, Ji] = [H, J2] = [H, A, ] = [H, B]=0): 

d d d 
J^ = -i{uo- ^1^), J2 = -i{u2^ ""3^), (4.4) 

OUi OUq OU^ OU2 

I f 2 o ^ ^ . "V Ci , C2 

A = — - u -—^ — UjUi- — 6Ur 



4 \ du^ duiduj dui J 2 \Uq + Ui M2 + % 

(4.5) 

_^ Ci C2 



We then have the polynomial algebra ( with J3 = —\Ji — \J2 and J3 = 
-i Ji + I J2 and [A, J3] = [A, 4] = [B, J3] = [S, J'^ = 0) 

[A, C] = 2{A, B} + AJ^f^H - 2(ci - Ca)//, (4.7) 
[B, C] = -2B^ - IGh^u^A + 2H^ - 8a;V| - SwV^^ ^ g^^^ ^ _ 1)^2^ 

for which the structure constants may be obtained by comparing (2.1) and 
(4.7) : 

(3 = a = S = e = 0, 7 = 2, d= -16a^ (4.8) 
C = U3J3H - 2(ci - C2)H, z = 2H^ - 8a;2(j| + Jg^) + 8(ci + ca - 1)uj^ 
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We can obtain the Casimir operator and rewrite it in terms of H, J3 , J3 (which 
commute with all other integrals of motion) : 

K = -AJlH^ - AJ!^H^ + 4(ci + C2)H^ + ^oJ^Jp'^^ (4.9) 

-16(ci - C2)a;V3 + 4(ci - C2) V. 

The structure constants (4.8) and the Casimir operator (4.9) allow us to 
obtain the structure function $(a;) of the deformed oscillator algebra. The 
energy and the two angular momenta are fixed, Hil) = Ei/j, J^ij) = mij) and 
JgV^ = sip. The structure function is then given by (with mi = (m + s)^ + 2ci 
and 1712 = {m — s)^ + 2C2) 

$(a;) = 3 ■ 222a;2 ( x + u -(--—]]( x + u -(- + —]] (4.10) 



2 2uj J J \ \2 2uj^ ^ 
X + u — [ {nil + "^12)]] ( X + u — ( (mi — m2' 



2 2' 7/ V V2 2 
X + u — ( {—mi + m2) ) ] ( X + u — ( (—mi — m2) 



2 2' 7/ V V2 2 

We found the following finite dimensional unitary representation : 

^ / mi + m2 , IE / , X 

E = 2u;{p+1+ u = --—, (4.11) 

4){x) = 3 ■ 2^^u^x{p + 1 - x){p + 1 + mi - x){p + 1 + m2 - x) (4.12) 

{p + 1 + nil + "^2 — x){2p + 2 + mi + m2 — x). 

In Ref.26 the systems was studied using double polar coordinates and the 
double polar quantum numbers ni and 77-2. There are relations similar as 
Eq.(3.13). These results coincide by taking p = ni + n2- 

To summarize : We have presented an algebraic derivation of the energy 
spectrum of this superintegrable system, and its finite- dimensional unitary 
representations. This system is interesting not only for the duality (4.1) but 
also because 4D systems have not been studied systematically, and the appli- 
cation to them of polynomial algebra, with its power to determine spectra, 
is an unexplored subject. We leave to a future article the study of the corres- 
ponding classical superintegrable system and the polynomial Poisson algebra. 
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5 MICZ-Kepler system on 5^, deformed so(4) 
algebras and quadratic algebra 



The quantum mechanical problem of motion of a arbitrary dyon field modi- 
fied by a centrifugal term proportional to the square of the magnetic charge, 
(the MICZ-Kepler problem) has also been considered in the 3D space of 
constant positive curvature [23]. The integrals of motion of this superin- 
tegrable system form a cubic algebra similar to that of the Kepler-Coulomb 
problem on S^. The explicit form of this cubic algebra may be used to obtain 
the energy spectrum of the problem [23] . We will point out that we can form 
a quadratic algebra with three generators and use results obtained in context 
of 2D systems to obtain the energy spectrum. 

We begin by recalling some known results concerning the MICZ-Kepler 
system on S^. In coordinates x = {xi, X2, X3} and x^x'^ = x."^ + x1 = R"^ , the 
metric of the sphere is ds'^ = {dxY + ^iid the Hamiltonian 



//„^-i(V" + ,^m + .^) + ^-|^. (5.1) 



^1 = ^ 1 \(\ I I V' ^2 = ^ 1 |.| I . — T, ^3 = 0. (5.2) 

From [22] we then have 

TTa = -ida + A^, ^4 = -^^4, = X^-Ra - XaPi, (5.3) 

with 



X 

[TTa, Xb] = -iSab, ['rra,P4] = 0, [tT^, TTfe] = i^eabc-r4:T, [P4, X4] = -^, (5.4) 

\xr 



[Na, Nb] = leabcJc + R'Fab, [Ja, Nb] = iCabcNc, (5.5) 
[Ja, Jb] = i^abcJc, Fab = [7^a,'^b]- 

The integrals of motion are the Poincare and Laplace-Runge-Lenz vector 

Ja = eabcXbT^c- K = ^(^abc{JbN c - NbJ c) + ■ (5.6) 

\x\ in \x\ 
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The MICZ-Kepler Hamiltonian on the sphere can be rewritten 

2i?2 ^2/22|x|2 R\xy 
and the cubic algebra [21] is then given by 

[A,, A] = -1i{E - — + ^)ea6c Jc, (5.8) 

[•-^aj^ft] = i^abcAc, [Ja, Jb] = i^abcJc 

with the relations 

= 2H{J^ - /i^ + 1) - J2(J2 _ ^2 ^ 2) + (5.9) 

1- J= -a//. (5.10) 



5.1 Quadratic algebra and realizations in term of defor- 
med oscillator algebras 

The representations of this algebra were discussed in terms of the cubic 
algebra [22] given by (5.8). Deformations of the Lie algebra so(4), so(3, 1) 
and e(3) that leave their so(3) subalgebra undeformed and preserve their 
coset structure were studied by Quesne [27] and the cubic algebra (5.8) is a 
particular case. We will point out how results obtained for polynomial algebra 
in the context of 2D superintegrable systems can be used to obtain the finite- 
dimensional unitary representation and the corresponding degenerate energy 
spectrum. The duality transformation for the Kepler-coulomb system in 2,3 
and 5 dimensional sphere was investigated in Ref. 28. Let us consider the 
three generators 

A = J\ B = A3, C=[A,B], (5.11) 
and form the quadratic algebra 

[A,B] = C, [A,C] = 2{A,B} + 4afiJ,, (5.12) 
[B, C] = -2B' - ^^A' + (SH + -1(J| + /.^ - 1))A 
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H2a' + AH{l-Jl-^^^)-2^-^). 
By comparing (2.4) and (5.12) we obtain the structure constants 

/3 = 5 = e = 0, 7 = 2, C = 4aM, « = (5-13) 

= 8i7 + -i( + /i^ - 1), z = 2a^ + AH{1 - - /i^) - 2^^. 

The structure constants can be used to obtain the Casimir operator K of 
the quadratic algebra. We rewrite this operator as a polynomial of the Ha- 
miltonian and using (2.2) and the two Casimir operators of the quadratic 
algebra (5.9) and (5.10) : 

t2 2 

K = Aa^fx^ + 4J>2 - 8-^ - SJ^.n^H. (5.14) 

The energy and the z-component of the Poincare vector are again fixed, 
H = E and = m. For the MICZ-Kepler systems on we then use (2.4), 
(5.13) and (5.14) to obtain the structure function 

$(x) = -^^{Am^-{l-2{x+u)f){-16R^a^-8ER\l-2{x+u))^+{l-2{x+u)f 

(5.15) 

+ {l-2{x + u)f{-?>-A{x + u)+A{x + uf)){l-A{x + u)+A{x + uf-Ai/). 
This can be rewritten in the following form (with E' = 1 + 2ER^) : 

3218 1111 
*(^) = — ^(^+^-(2(l-2^)))(^+^-(2(l+2m)))(z+«-(-(l-2Ai)))(2;+n-(-(l+2/x))) 

(x+u-{^-^^JE' - v^O' + AR'^a'')){x+u-{^+^^jE' - ^R^a^)) 

{x+u-{^-^^/e' + v^')' + ^R^a^)){x+u-{^+^^E' + ^R^a^))- 
We obtain, with the constraint given by (2.19), 
N"^ -1 1 

i? = -^ + ^^, u=-{l + 2\^i\), N = p+l + \fx\, (5.16) 
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2218 R^a^ 
^{x) = — ^x(p+l-x)(x+2|/i|)(x+m+|/i|)(x-m+|/i|)(p+l+x+|^|)((2x+2|^|)24 



(5.17) 

To summarize : We have obtained the finite dimensional unitary representa- 
tions and the corresponding energy spectrum, and have thus obtained a new 
algebraic derivation of the energy spectrum of this system using algebraic 
results obtained in the context of 2D superintegrable systems. 



6 Conclusion 

In this paper we considered the generalized MICZ-Kepler systems and the 
4D singular harmonic oscillator. These systems are respectively in the 3D 
and the 4D Euclidean space. We presented the quadratic algebra for these 
two systems (for which, to our knowledge, no polynomial algebras have been 
obtained before), their realizations in terms of deformed oscillator algebras, 
their finite dimensional unitary representations and energy spectrum. The 
use of polynomial algebras to solve superintegrable systems in dimensions 
higher than two is a relatively unexplored subject [29], and these systems, 
with the duality relation between them, provide an important example. 

We pointed out also that results obtained in [7-9] in the context of 2D 
superintegrability may be applied to other superintegrable systems in hi- 
gher dimensions. The structure constants of the quadratic algebras and the 
structure function involve not only the Hamiltonian and but also other in- 
volutive integrals. In the case of the generalized MICZ-Kepler system there 
is one other such integral of motion (J^) and in the case of the 4D singular 
oscillator there are two such integrals (Ji and J2). 

We also obtained a quadratic algebra for the MICZ-Kepler systems on 
from the cubic algebra with six generators and two Casimir operators. We 
used the results of Daskaloyannis [9] to obtain realizations of this in terms of 
deformed oscillator algebras and presented an alternative algebraic derivation 
of the energy spectrum. Thus the problem of finding the finite dimensional 
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unitary representations and the corresponding energy spectrum for this sys- 
tem can be achieved by applying the constraint (2.5) to the structure function 
and solving an algebraic system of equation. 

Similarly as for the Kepler- Coulomb systems, the Hartman potential and 
the Kepler- Coulomb system and harmonic oscillator on the three sphere 
[7,8] the the MICZ-Kepler systems on S^, the generalized MICZ-Kepler sys- 
tems and the 4D sungular oscillator possess a quadratic algebra with three 
generators and we can reobtain the energy spectrum from this algebra. Ho- 
wever, results obtained in context of two-dimensional systems should be ap- 
plied carefully in context of higher dimensional systems. A systematic study 
of higher dimensional suprintegrable systems and their polynomial algebras 
remains to be done. 

A systematic search for systems with spin was also started recently [30]. It 
could be interesting to see how results obtained in [7-9] could be applied in 
this context. 

The MICZ-Kpler systems [31-32], the MICZ-Kepler systems on the three 
sphere [33] and the generalized Kepler system [34] were discussed from the 
point of view of supersymmetric quantum mechanics. The factorization of the 
generalized MICZ-Kepler systems were also discussed [35]. We pointed out 
in context of two-dimensional superintegrable systems a relation between 
higher order integrals of motion, higher order ladder operators and higher 
order supersymmetric quantum mechanics [10,11]. Such relation should be 
to be more studied in context of systems with and without monopole. 
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